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ABSTRACT 

Experimental  investigations  were  undertaken  to  determine  the  collapse 
pressures  of  relatively  thick  circular  cylindrical  shells  subjected  to 
hydrostatic  pressure.  A  total  of  twenty- three  fiber  glass  reinforced 
plastic  cylinders  with  thickness  to  diameter  ratios  ranging  from  0.05  to 
0.09  was  tested  to  determine  collapse  pressures  and  the  effective  modulus 
of  elasticity  of  the  material.  The  length  to  diameter  ratio  was  kept  con- 
stant. Results  were  plotted  non-dimensionally  with  P/E  as  the  ordinate 
and  h/D  as  the  abscissa.  The  results  were  compared  with  available 
theoretical  instability  formulas  and  were  in  very  good  agreement  with 
theory  for  ratios  of  h/D  up  to  about  0.07.  Above  this  value  of  h/D, 
experimental  collapse  pressures  were  lower  than  those  predicted  by 
theory. 

It  was  concluded  that  the  von  Mises  equation  [equation  (6),  Ref.  22J 
is  the  best  instability  equation  for  cylinders  with  ratios  of  h/b  greater 
than  0.05.  It  was  also  concluded  that  for  almost  all  practical  cases, 
the  thin  shell  theory  is  sufficiently  accurate  for  predicting  instability 
failures. 

Derivations  for  the  following  are  included:  (1) Expressions  for 
calculation  of  the  moduli  of  elasticity  in  the  circumferential  and 
longitudinal  directions;   (2) An  expression  for  the  shear  modulus  and 
for  the  effective  modulus  of  elasticity  of  an  orthotropic  material; 
(3) A  determinant  expressing  the  critical  pressure  of  simply  supported, 
orthotropic,  circular  cylinders  subjected  to  hydrostatic  pressure;  and 
(A) An  expression  for  calculating  the  bending  stress  in  a  simply  supported 
circular  cylinder  subjected  to  hydrostatic  pressure. 
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NOTATION 
D     Flexural  rigidity,   D  = 


/2(/-z/2j 
D     Mean  diameter 

D^  Inside  diameter 

T>2  Outside  diameter  ** 

E  Modulus  of  elasticity 

Eq  Effective  modulus  of  elasticity  for  anisotropic  material 

Ex  Modulus  of  elasticity  in  the  longitudinal  direction 

Ee,  E^  Modulus  of  elasticity  in  the  circumferential  direction 

G  Shear  modulus,  G  =  zrr \ 

h  Thickness  of  shell 

Kj  Non-dimensional  term  from  thick  cylinder  stress  formula, 

1  =  Rl-R? 

K2    Non-dimensional  terra  from  thick  cylinder  stress  formula, 

K2  =-2£±5l  =  /<,-/ 
Rt  -  R,1 

L     Effective  length 

L-oa   Actual  length 

m     Number  of  waves  or  lobes  along  effective  length  of  cylinder  at 
time  of  collapse 

m^  Strain-pressure  ratio,  circumferential  strain  gage 

in*  Strain-pressure  ratio,  longitudinal  strain  gage. 

Mx  Axial  bending  moment  per  unit  of  circumference 

M^  Circumferential  bending  moment  per  unit  of  length 

MX0  Twisting  moment  per  unit  of  length  on  element  cut  out  by  generators 

M^x   Twisting  moment  per  unit  of  circumference  on  element  cut  out  by 
plane  perpendicular  to  axis 

## 

In  the  notation  used  in  this  paper,  the  subscript  "l"  refers  to  the 
inside  of  the  cylinder,  and  the  subscript  "2"  refers  to  the  outside  of 
the  cylinder. 


iv 


n  Number  of  lobes  or  waves  in  a  complete  circumferential  belt  around 
cylinder  at  time  of  collapse 

Nx  Normal  force  per  unit  of  circumference  in  axial  direction 

Np  Normal  force  per  unit  of  length  in  circumferential  direction 

Nx*  Shear  force  per  unit  length  on  elements  cut  by  generators 

N^x  Shear  force  per  unit  of  circumference  on  elements  cut  by  planes 
perpendicular  to  axis 

p  Pressure 

Qp  Radial  shear  in  the  xz  plane 

Qx  Radial  shear  in  the  yz  plane 

r  Radial  coordinate 

R  Mean  radius 

R^  Inside  radius 

R2  Outside  radius 

u  Axial  displacement 

v  Circumferential  displacement 

w  Radial  displacement 


z     Perpendicular  distance  from  mid-surface  of  shell,  positive  inward 

J-  Jl 
3'  Dx 


oc     Non-dimensional  terra  in  buckling  equations,  ot  = -^-  " 


T"  Shearing  strain 

Y"ij  Shearing  strain  along  i-plane  in  j-direction 

Cx  Normal  strain  in  axial  direction 

£0  Normal  strain  in  circumferential  direction 

€r  Normal  strain  in  radial  direction 

9  Cylindrical  (angular)  coordinate 

tt  D 

\  Non-dimensional  term  in  buckling  equations,      \=  m  ~2~L~ 

2/  Poisson's  ratio 


0TX  Normal  stress  in  axial  direction 

O^  Normal  stress  in  circumferential  direction 

7^  •  Shear  stress  along  i-plane  in  j-direction 

<P  Cylindrical  (angular)  coordinate 

*X  *  Change  of  curvature  of  generator 

%jt  Change  of  curvature  of  circumference 

%x0  Twist  of  element  due  to  Mxp 
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INTRODUCTION 

The  absence  of  experimental  results  and  elastic  instability  analyses 
for  thick  cylinders  has  been  pointed  out  by  Wenk-*-.  On  the  basis  of  this 
information  and  consultation  with  Mr.  John  Pulos  of  the  Structural  Mechanics 
Laboratory,  David  Taylor  Model  Basin,  the  authors  decided  to  undertake  an 
experimental  investigation  of  the  elastic  instability  of  relative] y  thick 
cylinders  subjected  to  external  hydrostatic  pressure.  Thus,  the  primary 
purpose  of  this  thesis  evolved;  which  is  the  experimental  evaluation  of  avail- 
able theoretical  instability  formulas  as  applied  to  cylinders  with  thickness 
to  diameter  ratios  greater  than  those  heretofore  tested.  From  the  results 
of  the  proposed  experiments  it  was  anticipated  that  the  following  question 
posed  by  Dr.  Wenk  (see  "Discussion"  of  Ref.  1)  could  be  answered: 

"At  what  ratio  of  thickness  to  diameter  are  errors  in  the  thin-shell 
approximation  unacceptable?" 

A  prime  requisite  prior  to  undertaking  the  thesis  was  the  making  of  a 
careful  survey  of  the  literature  to  determine  the  extent  of  theoretical  and 
experimental  work.  The  authors  conducted  such  a  survey  both  prior  to 
beginning  the  experiments  and  on  a  continuing  basis  throughout  the  period 
of  the  thesis  work.  The  results  of  this  search  are  presented  in  the  section 
entitled  "Literature  Survey." 

It  was  also  necessary  to  design  and  construct  the  experimental  appa- 
ratus and  the  models.  The  materials  and  methods  of  design,  together  with 
the  procedures  employed,  are  indicated  in  the  body  of  this  paper.  Since  the 
models  were  fabricated  from  glass -reinforced  plastic,  the  e]astic  formulas 
required  in  analyzing  the  data  were  modified  to  account  for  the  difference 
in  the  values  of  the  principal  moduli  of  elasticity.  Derivations  for  the 
expressions  used  are  presented  in  the  appendix. 

1 
References  are  listed  beginning  on  page  67. 
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The  thesis  begins  with  a  section  on  concepts  and  assumptions  employed 
in  the  theory  of  thin  shells.  This  study  was  made  in  order  to  gain  an 
insight  into  the  differences  between  a  thick  and  a  thin  shell  theory. 
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CONCEPTS  AND  ASSUMPTIONS  IN  THE  THEORY  OF  THIN  SHET.T.S 

A  thin  shell  may  be  defined  as  one  for  which  the  ratio  of  the  thick- 
ness of  the  shell  to  the  principal  radii  of  curvature  is  small  compared 
to  unity.  The  problem  of  classifying  a  shell  as  thick  or  thin  therefore 
reduces  to  the  determination  of  a  value  of  h/R  below  which  solution  errors 
will  be  within  a  prescribed  limit.  Novoxhilov  has  stated  that  for  or- 
dinary technical  calculations  it  is  admissable  to  have  a  relative  error 

of  5%,   and  that  for  thin  shells  this  corresponds  to  a  maximum  h/R  =  0.05. 

3 
Timoshenko  (p. 354)  has  suggested  a  value  of  h/R  =  0.10  as  the  dividing 

point  for  thin  curved  bars. 

There  are  three  principal  assumptions  employed  in  the  theory  of  thin 

shells  in  addition  to  those  made  for  most  problems  in  elasticity.  They  are 

(1)  The  normal  stresses  acting  on  planes  parallel  to  the  middle 
surface  may  be  neglected  in  comparison  with  the  other  stresses. 

(2)  The  normals  to  the  undeformed  middle  surface  are  deformed  with- 
out change  of  length  into  the  normals  to  the  deformed  middle 
surface. 

(3)  The  thiokness  h  is  very  small  compared  to  the  least  radii  of 
curvature  (i.e.  h/R«l). 

Assumptions  (1)  and  (2)  essentially  reduce  the  problem  to  the  case  of 

plane  stress.  Consider  the  generalized  matrix  form  of  Hooka's  law  (see 

Borg  ,  p.  59). 


where 


A  = 


T  -Xl^     +  2/url 


T  =  the  stress  tensor 
1   =  the  strain  tensor 
I]=  the  trace  of  the  strain  tensor 
Zy=   the  unit,  3x3,  matrix 
X^/a  =  constants 

^K _     E 
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Now  at  the  inside  surface  of  the  cylinder  (i.e.  at  z  =  h/2)  the 
stress  component  o~z  is  aero.  Then,  if  the  shell  ia  thin,  this  compo- 
nent can  be  assumed  to  be  zero  throughout  the  thickness  of  the  shell 
since  it  is  small  compared  to  the  stresses  <rm   and  <r_. 

Assumption  (2)  is  equivalent  to  assuming  that  the  shearing  strains 

V   and  YC  .  are  zero,  for  the  definition  of  shearing  strain  is  the 
xz     e  z 

change  in  value  of  an  originally  right  angle  from  the  unstrained  state. 
The  assumption  is  that  this  right  angle  is  preserved  during  deformation. 

Therefore,  letting  0~z  =  ^xz=^zDe  zero,  the  tensor  equation 
becomes 


0h    <Ke      OZz 


=  A(6,^e  +  €x) 


I     0     0\ 

0      I      0 


€Kf  M 


+   Z/u      fcVS 


i<r*z     <rez     o  I  \0     O     / J  \0 


e 


Expanding  the  equation  gives 

ai    =   O  =   A  (€*  +  €0  i-  €z  )  *■  Z^u  €x 

<r*z  =  <^2  =  o 

Eliminating  £  z  from  the  first  two  equations  by  use  of  the  third  gives 
the  desired  stress-strain  relations 

<Z   =  /f*,*     (€*   **^e) 


Furthermore,  assumptions  (1)  and  (3)  are  utilized  in  reducing  the 
strain  components  to  a  usable  form  (see  Wang'9,  pp.  335-34.0).  When  this 


-  4-- 


is  done,  the  equation  for  cr     becomes 


K    = 


/-  >"z 


fc*o    ^^eo-z(X, 


-h      V 


x*)] 


where  <)n 


xK   = 


<2. 


X0  -  a. 


^e  (  c^ 


J-   i+£ 
a,    ~f& 


) 


The  expressions  for  the  other  stress  components  are  similar  in  degree 
of  complexity. 

Finally,  to  obtain  the  relations  between  the  stress  resultants 
and  couples,  the  following  expressions  must  be  integrated. 


Ax 


('+%*)* 


Again,  assumption  (3)  is  used  and  z/R»  is  neglected.  Then,  Nxe  =  Nex, 
Mjq  =  Mgjj,  and  the  expressions  may  be  integrated  without  difficulty. 
For  example,  the  expression  for  Nx  becomes 


This  equation  and  the  remaining  ones  obtained  in  a  similar  manner  are 
known  as  Love's  first  approximation. 

It  is  possible  to  derive  expressions  for  the  stress  resultants  and 

couples  without  resorting  to  the  use  of  assumptions  (1)  -  (3).  Naghdi  , 

using  a  variational  method  due  to  Reissner  ,  has  made  such  a  derivation. 

In  order  to  gain  an  appreciation  for  the  simplicity  of  the  above  expression 

and  of  the  complexities  involved  in  a  thick  shell  theory,  the  more  exact 

expression  for  Nx  is  given  below. 
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Presumably,  these  more  exact  relations  could  be  Introduced  into  the 
equilibrium  equations,  and,  provided  a  solution  could  be  effected,  a 
thick  shell  analysis  for  the  instability  problem  would  become  available. 
An  analytical  approach  to  a  thick  shell  instability  analysis  has  not 
been  attempted  in  this  thesis,  however. 
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LITERATURE  SURVEY 

The  theoretical  and  experimental  analysis  of  the  buckling  of 
cylindrical  shells  subjected  to  hydrostatic  pressure  has  attracted  the 
attention  of  numerous  investigators  dating  from  experiments  conducted  by 
Fairbairn  in  1858.  All  investigators  employed  the  theory  of  thin  shells 
in  arriving  at  theoretical  formulas  and  most  experiments  were  carried  out 
on  models  with  relatively  thin  walls.  Tables  and  descriptive  information 
are  presented  in  this  section  to  indicate  the  extent  of  previous  theoretical 
and  experimental  investigations.  Particularly  helpful  in  the  making  of 

Q 

the  literature  survey  were  two  reports  by  Nash7  which  contain  a  listing 
of  2339  papers  and  books  relevant  to  experimental  and  theoretical  work 
on  shells  and  shell-like  structures  from  1828  through  1956. 


Results  of  Previous  Experimental  Investigations 

Table  1  is  a  chronological  listing  of  various  investigators  together 
with  the  range  of  geometries  involved  and  other  pertinent  data.  Since  the 
first  analytical  solution  to  the  instability  problem  was  not  available 
until  1913,  all  investigations  prior  to  this  time  were  made  principally 
to  determine  empirical  design  formulas.  Those  investigations  made  in  later 
years  were  carried  out  to  test  the  validity  of  analytical  solutions  either 
in  general  or  as  applied  to  a  particular  problem. 

With  the  exception  of  Southwell,  all  experimental  investigators  prior 
to  1929  restricted  the  loading  on  the  cylinder  to  a  radial  pressure  only. 
The  axial  thrust  was  eliminated  by  the  use  of  a  relatively  stiff  rod  which 
connected  the  end  plugs.  The  use  of  radial  pressure  only  resulted  from  the 
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absence  of  analytical  formulas  based  on  combined  loading,  and  from  the 
fact  that  one  of  the  principal  applications  of  the  results  was  to  the 
study  of  boiler  flues.  Southwell  did  use  a  combined  loading  to  simplify 
the  apparatus.  He  neglected  the  effect  of  axial  loading,  however,  in 
comparing  test  results  with  his  analytical  approach. 

Most  investigators  after  1929  employed  models  subjected  to  hydrostatic 
(combined  radial  and  axial)  loads.  An  analytical  solution  to  the  problem 
of  the  cylinder  under  hydrostatic  pressure  was  available  by  this  time,  and 
the  submarine  presented  an  immediate  application. 

From  columns  4,5,  and  6  of  Table  1  the  following  is  noted:  Most  of  the 
models  were  relatively  small  in  diameter,  as  was  the  h/D  ratio.   In  those 
cases  where  the  h/D  ratio  was  greater  than  0.05  the  model  was  either  so 
long  that  the  values  of  the  critical  pressure  approached  that  of  an  infinite 
tube,  or  they  were  in  the  range  where  yielding  of  the  material  occurred 
rather  than  an  instability  failure. 

The  type  of  support  used  hy  the  various  investigators  is  indicated 
in  column  7.   In  general,  this  information  was  procured  from  sketches  of 
the  models  as  prepared  for  testing.  A  clamped  type  of  support  is  indicated 
only  where  the  author  concerned  specifically  labelled  the  test  condition  as 
such.  It  is  again  noteworthy  that  no  attempt  was  made  to  clamp  the  ends 
until  1941   when  an  analytical  solution  became  available.  Also,  as  indi- 
cated by  Cook",  simply  supported  ends  can  be  approached  experimentally 
much  more  closely  than  fixed  ends. 

Column  8  indicates  the  maximum  collapse  pressures  achieved  by  the 
various  investigators;  column  9  the  number  of  models  tested;  and  column 
10  the  variety  of  materials  which  have  been  utilized  in  the  construction 
of  models. 

Several  general  observations  can  be  made  based  on  the  survey  of 
experimental  results  found  in  the  literature: 
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(1)  Scale  effect  is  not  important  provided  reasonable  care  is  taken  in 
fabricating  the  models.  After  testing  100  models  representing 
the  strength  hull  of  a  submarine,  Windenburg^  concluded  that 
scale  effect  was  not  evident.  Some  the  earlier  experiments 

have  not  given  good  results  when  compared  to  theory  because 
of  the  inadequacy  of  materials  and  unsound  fabrication  techniques. 
Most  of  the  results  obtained  in  the  last  two  decades  have  given 
very  good  correlation. 

(2)  Although  the  assumed  end  condition  is  taken  as  either  simply 

supported  or  as  clamped,  neither  condition  can  be  precisely  ob- 

13 

tained  in  practice.  Based  on  tests  of  258  brass  models,  Cornell  J 

concluded  that  the  end  plugs  inserted  into  the  ends  of  the  tubes 
offered  a  restraint  that  tended  to  reduce  the  effective  length. 
He  found  that  to  get  close  agreement  of  his  test  results  with 
theory,  it  was  necessary  to  use  an  effective  free  length  which 
differed  from  the  actual  free  length  by  ^%   to  1%,      (Because  of 
the  small  value  of  such  a  correction,  the  authors  assume  the 
measured  free  length  to  be  the  same  as  the  effective  free  length 
for  the  tests  reported  in  this  thesis.) 

(3)  Several  investigators  have  conducted  carefully  controlled  experi- 
ments to  determine  the  influence  of  other  factors,  such  as  eccen- 
tricity, on  the  collapse  pressure.  Sturm^  was  probably  the  first 
to  investigate  the  problem  of  out-of-roundness  from  an  experimental 
and  analytical  point  of  view.   In  recent  years,  Cleaver^-  has 
conducted  very  thorough,  extensive,  and  carefully  controlled 
tests  on  530  models  subjected  to  radial  pressure.  He  lists  the 
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following  factors  as  affecting  the  collapse  pressure: 

a.  Eccentricity  of  bore  relative  to  the  external  surface. 

b.  Variations  from  true  circular  shape  (out-of-roundness) . 

c.  Ovality  of  the  specimen  (the  difference  in  the  maximum  and 
minimum  diameters  divided  by  the  nominal  diameter) . 

d.  Variation  of  material  properties,, 

From  the  results  of  his  tests  Sturm  concludes  that  eccentricity  of  bore 
is  the  predominant  manufacturing  variable  affecting  collapse  strength, 
but  for  eccentricities  within  the  limits  imposed  by  current  specifications, 
its  effects  are  small  and  can  be  neglected  for  practical  purposes. 
Variations  in  circularity  exert  no  measurable  systematic  influence,  and 
provided  they  are  within  current  specifications  (i.e.  for  commercial  tubes), 
they  too  may  be  neglected. 

The  authors  point  out  in  passing  that  for  a  well-made,  simply 
supported  model  the  above  effects  tend  to  cancel  the  effect  of  any 
end  restraint  that  might  be  present.  Therefore,  neglecting  these  effects 
may  be  justified  when  attempting  to  correlate  experimental  and  theoretical 
results. 


Results  of  Analytical  Investigations 

Table  2  presents  a  summary  of  results  of  principal  theoretical 
investigations.  The  method  of  approach  to  the  problem  is  indicated. 

The  differences  in  the  various  formulas  can  be  accounted  for  by 
the  inclusion  or  omission  of  certain  higher  order  terms  in  the  differ- 
ential equations  of  equilibrium,  by  the  method  of  approach  (i.e.  equil- 
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ibriura,  energy,  kinetic  ),  and  by  the  degree  of  simplification  of  the 

final  result. 

All  of  the  more  exact  formulas  are  found  to  have  the  form 

p/B  =  Fx   (h/D)3  +  F2  (h/D) 

where  F.  and  F„  are  functions  of 
1     2 

(a)  The  length  to  diameter  ratio. 

(b)  Poisson's  ratio. 

(c)  The  number  of  lobes  in  a  complete  cicumferential  belt  at  the 
time  of  collapse. 

(d)  The  number  of  lobes  along  a  generator  of  the  cylinder  at  the 
time  of  collapse.  Von  Mises1*  has  shown  this  quantity  to  be 
unity  for  minimum  buckling  pressure. 

In  order  to  facilitate  comparison  of  the  buckling  equations  with 
test  results,  values  of  F.  and  F2  for  the  model  geometry  and  material 
under  consideration  in  this  paper  have  been  calculated  for  the  various 
buckling  formulas.  These  values  are  given  in  Table  3.  Curves  for  p/E 
versus  h/D  are  also  included  for  these  cases.  (See  Fig.  13). 

In  addition  to  the  formulas  for  the  buckling  pressure  of  cylinders 
under  hydrostatic  load,  it  is  of  interest  to  examine  formulas  for  buck- 
ling pressure  under  uniform  radial  load  only.  Table  U   gives  some  of  the 
principal  radial  load  buckling  equations.  A  comparison  of  the  radial 
load  formulas  with  the  hydrostatic  load  counterparts  indicates  that  the 

following  relationship  is  approximately  true 

2 

buckling  pressure  under  hydrostatic  load  =    n  -  1 


buckling  pressure  under  radial  load        tr   -  1  +  A 

2 

Thus  as  A   approaches  aero  (L/D  -**),  the  formulas  for  hydrostatic 

load  should  yield  the  same  results  as  the  formulas  for  radial  load. 

For  the  case  where  A    -   0,  n  becomes  two,  and  all  of  the  more  exact 

hydrostatic  and  radial  load  buckling  pressure  equations  reduce  to  the 

Bresse  equation  (Table  4-)  for  a  tube  of  infinite  length. 
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TEST  APPARATUS.  MODELS,  AND  TEST  PROCEDURE 

Test  Apparatus 
The  test  apparatus  consists  of  a  pressure  vessel,  pump,  tubing, 
pressure  gages,  and  equipment  required  for  the  taking  of  strain  gage  data. 
The  arrangement  and  details  of  the  apparatus  are  shown  in  Figures  1  through 
5.  Specifications  which  were  drawn  up  for  the  test  apparatus  are  given 
below. 

Size  of  chamber:  5"  I.D.,  13"  depth 
Working  pressure:  10,000  psi 

Penetrations:  One  fitting  for  pressure  tubing,  one  fitting  for 
four  wires  from  strain  gages  on  the  external 
surface  of  the  model,  one  fitting  for  1/4."  tube 
to  vent  model  to  the  atmosphere  and  to  lead  out 
internal  strain  gage  wires. 
Pump:  Hand  operated  pump  with  self-contained  check  valve  and 
reservoir— the  type  commonly  used  with  hydraulic  jacks. 
Pressure  gages:  One  $000  psi  8W  Bourdon  gage  graduated  in  in- 
crements of  50  psi,  and  one  15,000  psi  6W  Bourdon 
gage,  graduated  in  increments  of  100  psi. 
Strain  Indicator:  Baldwin,  Type  N 
The  pressure  vessel  was  constructed  in  accordance  with  a  basic  design 
drawn  up  by  Mr.  Kenneth  Horn  of  the  David  Taylor  Model  Basin.  The  design 
was  modified  slightly  by  the1  authors,  to  suit  the  purpose  required.  The 
following  calculations  show  the  factors  of  safety,  for  a  working  pressure 
of  10,000  psi.  The  material  is  low  carbon  steel,  the  yield  strength  of 
which  was  assumed  to  be  30,000  psi.  Formulas  used  are  given  in  reference 
(39). 


-/?- 


t  =  wall  thickness,  inches,  =  f2'25  for  walls   ) 

(.2.125  for  bottom) 

d  =  inside  diameter  of  chamber,  inches,  =  5.00 

P  =  chamber  pressure,  psi 

S  =  maximum  stress,  psi,  =  30,000 

R  =  inside  radius  of  chamber,  inches,  =  2.50 

F.S.  =  factor  of  safety  based  on  yield  strength  of  30,000  psi 

Based  on  circumferential  stress  in  wall: 

P  *   st  *    =  (?0i0OO)(2.2?) =17,532  psi 

R  +  0.6t    2.50  +  (0.6)(2.25) 

F.S.  =  IZiil2  =  1.75 
10,000 

Based  on  longitudinal  stress  in  wall: 

P  =   a*    =  (2)(^»000).(2»^)n  =  8^,375  psi 
r   R  -  0.U     2.50  -  (0.4)(2.25)      '    P 

F.S.  =fH^  =  8-U 
10,000 

Based  on  stress  in  bottom: 

p  -         St2  _     (30.000)  (2. 125)2     _     -,,  j,Q       . 

P  =  d^(0.1o2)     "     (5.00)^0.1625  "     33»U9psi 

F  s  =  33,449  =  3.3z 
r,b*   10,000    ■'•■* 


*  Ref.  (39),  para.  UG-27(c)(l) 

*•  Ref.  (39),  Section  VIII,  para.  UG-34(1). 
Para.  UG-36(3)  states: 

"Single  openings  in  vessels  not  subject  to  rapid  fluctuations  in 
pressure  do  not  require  reinforcement  other  than  that  inherent  in  the 
construction  under  the  following  conditions: 

..».(b)  Threaded,  studded,  or  expanded  connections  in  which  the  hole 
cut  in  the  shell  or  head  is  not  greater  than  two-inch  pipe  size." 
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The  pressure  vessel  was  proof  tested  to  11,000  psi  on  30  August, 
I960.   It  was  raised  to  this  pressure  three  tiroes,  and  strain  data  were 
taken  on  the  outside  surface  using  a  circumferential  strain  gage.  This 
test  was  made  prior  to  the  making  of  the  penetrations  for  fittings  "A" 
and  "Bw  (Fig.  U) •  The  strain  data  taken  plotted  linearly  from  zero  to 
11,000  psi. 

This  apparatus  should  not  be  used  with  a  working  pressure  in  excess 
of  10,000  psi,  unless  the  pressure  vessel  is  again  proof  tested  at  the 
higher  pressure  with  strain  gage  measurements  taken  in  regions  of  high 
stress.  In  any  case,  it  is  recommended  that  a  pressure  of  15,000  psi 
not  be  exceeded,  and  that  for  working  pressures  in  excess  of  10,000  psi 
the  pressure  vessel  be  placed  in  a  pit  or  behind  a  suitable  barrier. 
All  tubing  and  fittings  used  should  be  rated  at  the  maximum  working 
pressure  employed. 

Persons  making  further  use  of  this  apparatus  are  warned  that 
valve  "G* ,   Fig.  4»  i£  rated  at  3000  psi.  It  is  recommended  that  this 
valve  be  replaced  by  one  rated  at  10,000  psi  or  higher.  As  a  safety 
precaution  in  performing  the  tests  for  this  thesis,  a  1/8"  steel  plate 
was  attached  to  the  table  between  the  operators  ajid  fittings  "D"  through 
"H"  (Fig.  4).  The  plate  can  be  seen  in  the  photograph,  Fig.  5. 
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FIGURE  5 
Photograph  of  Assembled  Apparatus 
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Model  Design 
Design  of  the  models  was  based  on  the  following  considerations: 

(1)  Maximum  working  pressure  of  test  apparatus:  10,000  psi 

(2)  Dimensions  of  test  chamber:   5"  I.D.  and  13"  depth 

Since  the  purpose  of  the  thesis  is  to  evaluate  the  instability  formulas 
for  relatively  thick-walled  cylinders,  it  was  necessary  to  design  the  models 
to  fail  by  instability,  and  to  have  a  relatively  large  thickness-diameter 
ratio.  Mr.  John  Pulos  of  the  David  Taylor  Model  Basin  recommended  a  range 
of  h/D  from  0.05  to  0.10.  Also,  Timoshenko^  (p. 354)  has  given  a  value  of 
h/D  =  0.05  as  an  estimated  division  point  between  thick  and  thin  curved 
bars.  Accordingly,  the  range  of  values  for  h/D  was  selected  as  0.05  to 
0.10. 

A  calculation  revealed  that  a  material  with  a  relatively  low  modulus 
of  elasticity  was  required  in  order  to  cause  instability  failures  in  the 
range  of  h/D  greater  than  0.05,  and,  at  the  same  time,  remain  within  the 
upper  limit  of  test  pressure  of  10,000  psi.  Captain  Hinners  suggested  the 
possible  use  of  glass-reinforced  plastic,  especially  in  view  of  the  current 
interest  in  this  material  because  of  its  favorable  strength-weight  ratio. 

The  David  Taylor  Model  Basin  had  published  a  report-*8  (No.  1413)  on 
the  hydrostatic  pressure  tests  of  a  cylindrical  shell  of  a  glass  fiber 
reinforced  epoxy  resin  fabricated  by  the  Zenith  Plastics  Co.  of  Gardena, 
California.  The  circumferential  modulus  of  this  material  was  reported  to 
be  4.8  x  10  psi,  and  the  plot  of  pressure  vs.  strain  was  linear  with  but 
a  slight  departure  from  linearity  near  the  collapse  pressure.  This  material 
seemed  to  meet  the  requirements  well.   It  had  disadvantages,  however,  in 
that  it  was  neither  homogeneous  nor  isotropic,  thus  violating  two  basic 
assumptions  in  the  usual  approach  to  the  theory  of  elasticity. 
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It  was  felt,  nevertheless,  that  these  disadvantages  would  be  overcome 
by  computing  an  "effective  modulus  of  elasticity"  from  strain  data,  and 
by  using  models  fabricated  in  such  a  way  as  to  give  maximum  dispersion 
of  the  glass  in  the  resin.  Correspondence  with  the  Zenith  Plastics  Co. 
was  initiated  to  obtain  information  regarding  the  procurement  of  the 
nedessary  models.  Zenith  very  generously  offered  to  supply  the  models  at 
no  cost. 

At  this  point  specifications  were  drawn  up  for  the  model  geometries. 
Using  DTMB  formula  (9)  [Formula  No.  5,  Table  2]  ,  a  set  of  curves  was 
drawn  for  various  values  of  L/D,  with  buckling  pressure  as  the  ordinate 
and  h/D  as  the  abscissa.  These  curves  are  shown  in  Fig.  6.  The  hoop 
stress  failure  pressure  was  also  plotted  on  these  curves  usinp  a  nominal 

value  of  compressive  strength  of  70,000  psi  as  obtained  from  DTMB  Report 

38 
LQ.3.    These  curves  were  then  used  to  select  appropriate  values  for 

L/D  and  h/D.   It  was  decided  to  keep  the  value  of  L/D  constant  at  about 

2.5  and  vary  the  value  of  h/D  from  about  0.05  to  0.09.  It  was  also 

decided  to  test  a  minimum  of  three  models  of  each  geometry  in  order  to 

determine  that  test  results  could  be  reproduced.  Table  5  lists  the  models 

received  from  Zenith  and  gives  their  measured  average  dimensions.  Fig.  7 

is  a  sketch  of  the  model  showing  the  end  closure  plugs.  The  plugs,  made 

from  medium  steel  stock,  were  designed  to  give  simple  support  of  the  model 

ends,  the  condition  for  the  von  Mises  buckling  equation.  The  effective 

length,  L,  is  as  shown.  The  inside  diameter  was  chosen  as  3.395"  because 

of  the  availability  at  the  Zenith  plant  of  a  mandrel  of  that  diameter 0 
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Model  Fabrication 
The  models  were  built  up  using  layers  of  glass  fiber  reinforced 
epoxy  resin  tape.  The  glass  fibers  ran  in  one  direction,  along  the 
length  of  the  tape  strip,  and  the  tape  was  applied  to  the  mandrel  both 
by  circumferential  winding  under  tension  (circumferential  plies)  and 
by  laying  on  strips  parallel  to  the  axis  of  the  mandrel  (longitudinal 
plies).  The  process  was  started  and  ended  with  two  circumferential  plies 
and  was  continued  using  alternate  layers  in  a  ratio  of  two  circumferential 
plies  to  one  longitudinal  ply.  This  2s 1  ratio  of  circumferential  to  long- 
itudinal plies  results  in  the  allowable  stresses  in  the  longitudinal 
and  circumferential  directions  being  in  close  agreement  with  the  induced 
stresses  in  these  directions.  The  resulting  structure  is  of  course  aniso- 
tropic. The  thickness  was  varied  by  varying  the  number  of  layers  of  tape. 
Five  tubes  were  fabricated,  each  about  60"  in  length  and  each  of  a  different 
thickness  (see  Table  5).  Upon  completion  of  the  curing  process,  most  of 
the  models  were  cut  from  the  tube  to  proper  length  while  still  on  the 
mandrel  at  the  plant.  Some  of  the  models  were  cut  from  the  remaining 
portions  of  tubes  after  delivery,  using  a  band-saw,  and  grinding  or 
machining  the  ends  until  smooth  and  normal  to  the  cylinder  axis.  By 
weighing  and  measuring  some  of  the  models,  the  density  of  the  material 
was  determined  to  be  about  0.07  pounds/cubic  inch.  Values  of  resin 
content  for  the  model  material  were  provided  by  Zenith,  and  are  given 
in  Table  11. 
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Instrumentation  of  Models 

Strain  data  were  taken  to  provide  a  means  for  calculating  an 
effective  modulus  of  elasticity  for  the  models.  As  shown  in  Appendix  C, 
values  of  the  strain-pressure  ratios,  taken  in  the  directions  of  principal 
stresses,  can  be  used  in  conjunction  with  Hooke's  Law  and  thick  shell 
theory  to  solve  for  E^  and  Ex.  A  value  for  Poisson's  ratio  must  be 
known  or  assumed. 

In  order  to  determine  average  values  for  e^  and  £„,  strain  gages 
were  applied  to  a  number  of  the  models,  oriented  both  circumferentially 
and  longitudinally.  The  number  of  gages  used  was  limited  both  by 
financial  considerations  and  by  the  number  of  leads  that  could  be  fed 
through  the  pressure  fittings  in  the  tank.  In  general,  four  gages  were 
applied  to  each  instrumented  nodel.  These  gages  were  located  at  the  mid- 
length  of  the  model  to  avoid  end  effects,  and  were  oriented  two  in  a 
circumferential  direction  and  two  in  a  longitudinal  direction.  The  two 
circumferential  gages  were  placed  "back  to  back",  inside  and  outside,  as 
were  the  two  longitudinal  gages. 

The  choice  of  gage  types  was  governed  by: 

(1)  the  need  to  have  non-pressure-sensitive  gages; 

(2)  the  necessity  for  measuring  strains  up  to  15,000  microinches/inch; 

(3)  the  desire  to  use  a  gage  length  long  enough  to  reduce  the  effect 
of  local  variations  in  strain;  and 

(4)  the  cost  of  gages. 

On  the  basis  of  these  considerations,  two  types  of  gages  were  selected. 
For  use  on  the  inside  of  the  model,  where  the  gage  would  not  be  subject  to 
the  pressure  of  the  fluid,  the  SR-A,  A-5-1  wire  gage  was  used.  This  gage 
has  a  1/2"  gage  length,  is  relatively  inexpensive,  and  is  suitable  for  the 
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measurement  of  strains  up  to  l£  to  2%,     For  use  on  the  outside  of  the 
model,  where  the  gage  would  be  subject  to  external  pressure  of  the  fluid, 
the  SR-4.,  paper-backed  constantan  foil  gage  FAP-50-12  was  selected.  This 
gage,  with  a  l/2n  gage  length,  is  not  as  pressure-sensitive  as  the  A-5-1 
gage.  It  also  is  suitable  for  measurements  of  strain  up  to  l£  to  2%, 

Compensating  gages  of  the  same  type  and  lot  number  were  used  in  con- 
junction with  the  active  gages. 

The  gages  were  applied  using  Duco  cement,  and  the  gage  manufacturer's 
instructions  were  followed.  After  the  cement  had  dried,  all  gages  were 
connected  in  the  strain  indicator  circuit  and  pressed  with  a  pencil  eraser 
tip  while  watching  the  strain  indicator  to  check  for  excessive  deflections 
of  the  needle. 
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Teat  Procedure 

Strain  gages  were  applied  to  the  inside  of  the  rodel  by  using  a 
"harness11— a  sheet  of  graph  paper  to  which  the  gages  were  attached.   This 
sheet  was  trimmed  to  exact  dimensions,  rolled  into  a  cylinder,  and  positioned 
inside  the  model  with  glue  applied  to  the  gages.  Leads  of  #28  stranded 
wire  were  attached  to  the  strain  gages. 

The  ends  of  the  model  were  prepared  by  thorough  cleaning  with  acetone. 
PR-1321  Class  A  sealant  *  was  then  applied  to  the  ends,  and  the  end  plugs 
inserted,  leading  the  wires  from  the  inside  of  the  model  out  through  the 
1/4"  tube  **  screwed  into  the  lower  end  plug.  After  the  sealant  set  up, 
the  model  was  ready  for  installation  in  the  tank.  Fig.  8  shows  a  model 
ready  for  installation  and  Fig.  5  shows  a  model  in  the  tank  and  partially 
lowered  into  testing  position.  Fig.  3  shows  a  model' inside  the  chamber. 

The  model  was  positioned  in  the  open  tank  by  first  leading  the  wires 
from  the  1/4."  end  plug  tube  through  the  "A"  fitting  in  the  center  of  the 
tank  bottom,  then  connecting  the  external  strain  gage  leads  to  the  leads 
inside  the  tank  from  the  "B"  fitting.  The  model  was  then  lowered,  allowing 
the  1/4."  tube  to  slide  through  the  "A"  fitting  until  the  lower  end  plug 
rested  on  a  steel  "spool"  support  about  2"  high  located  in  the  bottom  of 
the  chamber.  The  lock  nut  of  the  "A"  fitting  was  then  tightened,  taking 
care  not  to  over-tighten. 


*  Manufactured  by  the  Products  Research  Co.,  410  Jersey  Ave.,  Gloucester 
City,  New  Jersey 

**  This  tube  is  a  semi -permanent  attachment  to  the  end  plug,  its  threaded 
end  having  been  coated  with  sealant  prior  to  screwing  into  the  tapped 
end  plug. 

J 

•«*  It  was  observed  that  excessive  tightening  of  the  lock  nut  permanently 
deformed  the  ferrule  and  tube,  making  it  impossible  to  slide  the  tube 
back  out  of  the  fitting  after  the  test. 
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With  the  model  in  place,  the  strain  gage  leads  coming  from  the  l/Ln 
tube  and  from  the  "BH  fitting  underneath  the  pressure  vessel  were  connected 
to  the  strain  indicator  or  switching  unit.  The  chamber  was  filled  with  trans- 
former  oil  *  to  a  level  slightly  above  the  sealing  surface.  The  circular 
plate  was  inserted,  causing  the  excess  oil  to  be  squeezed  into  the  plug 
recess.  This  precaution  was  taken  to  insure  that  the  oil  was  relatively 
air-free.  The  screw  plug  was  then  positioned  and  tightened.  Tapping  the 
screw  plug  wrench  handle  with  a  hammer  insured  firm  seating  of  the  circular 

plate. 

Pressure  was  raised  by  operating  the  hand  pump.  The  valve  MGn  in  the 
line  to  the  5000  psi  gage  was  shut  at  4500  psi.  This  valve,  although  rated 
at  3000  psi,  performed  satisfactorily. 

In  some  cases,  only  one  strain  gage  was  measured  per  run,  using  the 
strain  indicator  directly.  In  most  of  these  tests,  at  least  two  runs  were 
made  per  gage,  and  the  readings  for  the  two  runs  were  found  to  vary  very 
little.  For  other  tests,  the  switching  unit  was  employed  and  a  number  of 
gage 8  were  read  per  run. 

Upon  completion  of  a  test,  the  screw  plug  was  removed,  and  the  "A" 
fitting  lock  nut  was  loosened  to  permit  removal  of  the  model.  As  the  model 
was  removed,  the  1/4"  tube  slipped  out  of  the  "A"  fitting,  allowing  the  oil 
in  the  tank  to  drain  into  a  receptacle  below. 

The  two  pressure  gages  used  during  the  tests  were  calibrated  at  the 
Material  Laboratory  before  beginning  the  tests  and  at  DTMB  toward  the  end 
of  testing.  Both  calibrations  were  made  using  dead-weight  gage 


*  WEMC0  Class  C  transformer  oil  was  U3ed  as  the  pressurizing  fluid,  making 
it  unnecessary  to  waterproof  strain  gages  exposed  to  the  fluid. 
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testers.  No  differences  were  noted  In  the  gage  calibrations.  The 
5000  psi  gage  was  adjusted  after  the  second  calibration  to  give  gage 
readings  closer  to  the  actual  pressure.  The  gage  calibration  curves 
are  shown  in  Figures  9  and  10. 
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EXPERIMENTAL  DATA 

A  total  of  twenty- three  models  was  tested,  twenty-one  of  which  were 
tested  to  failure.  Seven  of  the  models  were  instrumented  to  measure 
strains  in  order  to  determine  the  moduli  of  elasticity.  Three  models 
were  designed  to  fail  by  yielding  of  the  material.  This  was  done  in  order 
to  calculate  the  maximum  compressive  stress  which  the  material  would 
sustain  and  to  determine  the  appearance  of  models  after  yield  failure. 
Since  the  calculated  circumferential  stress  for  the  "D"  and  "E"  models 
was  close  to  the  nominal  strength  of  the  material,  this  information  was 
used  to  ascertain  whether  these  models  did,  in  fact,  fail  by  instability 
alone,  or  whether  some  yielding  might  have  taken  place. 

A  summary  of  the  models  tested,  with  the  respective  collapse  pressures, 
is  presented  in  Table  6.  The  recorded  strain  data  are  in  Appendix  A,  and 
the  pressure-strain  plots  are  in  Appendix  B.  Table  7  summarizes  the 
measured  strain-pressure  ratios  as  obtained  from  the  plots.  The  first 
models  instrumented  were  tested  with  the  compensating  gage  located  inside  the 
model  for  inside  active  gages,  and  with  the  compensating  gage  located  in  the 
pressurized  fluid  for  external  active  gages.  As  shown  in  Appendix  C,  it  is 
theoretically  possible  to  cancel  the  effect  of  radial  strain  on  the  outside 
active  gages  by  having  the  compensating  gage  mounted  on  a  block  which  is 
subject  to  the  hydrostatic  pressure.  Furthermore,  any  effect  of  pressure  on 
the  active  gage  in  the  fluid  would  be  compensated  for  ,  as  well  as  more 
accurate  temperature  compensation  provided.  The  later  models  tested,  however, 
had  all  gages — inside  and  outside — balanced  against  a  single  dummy  gage  in- 
side the  model.  This  latter  procedure  was  followed  because  it  had  been 
determined  that  the  error  caused  by  neglecting  the  Poisson  effect  of  the 
radial  strain  on  the  outside  gages  was  negligible,  and  therefore  the 
additional  complication  to  the  testing  procedure  was  unjustified. 
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In  order  to  be  able  to  compare  and  average  the  strain  data  taken  under 
the  two  different  methods  of  compensation,  the  strain-pressure  ratios  of 
active  gages  read  with  the  dummy  under  pressure  were  corrected  to  values 
equivalent  to  those  taken  with  the  dummy  gage  inside  the  model.  This  was 
accomplished  by  simply  adding  the  value  of  strain-pressure  ratio  of  a  gage 
mounted  on  a  dummy  block  under  pressure,  balanced  against  a  dummy  gage  inside 

a  model.  All  the  strain-pressure  ratios  used  in  computations,  therefore, 
have  values  which  would  result  from  active  gages  being  read  against  a 
dummy  located  on  the  inside  of  the  model,  not  under  pressure. 


--*/- 
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LAME  STRESS    IN  MODELS    AT    COLLAPSE    PRESSURE 
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A-2 

3/8o 

34,  9  lo 

17,  4SS- 

3/ ',    7  3o 

A-3 

31  30 

34,   16/ 

17,  /do 

31,   2.3/ 
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TABLE.       7 


MEASURED      STRAIN-  PRESSURE      (RATIOS 


MODE^ 

x   to*      [    %  to* 

rn  p           m*2 
X   \Ob          x  to* 

AV6-. 

™*XX  /O6 

AV<? 

A-  A- 

l_         . 

\    1.G2 

0.66,7 

/.  76       j  0.6,19 

A  -5              1/78      \    0  155     !    O.OUZ 
1 

(l-Z46,) 

(o.82  3)  !   18  4- 

A      AW& 

1-43 

| 
0-84S    \    ISO 

0.6/9 

1    I.6O 

1 

14  3 

O.I  3  Z 

C  -/ 

t.o6i    ;    0.76>o 

0.06,3      I./3I          o.  818         /.4&9 

0.798 

1.469 

1.131      \    0.8/3 

I 

\                 !                 1 

*> 

j 

i      D-2          \o.993 

j 

O  <c2S 

0.06,8    \(/.ob/)      (o.(o?j)>  I-Z38 

i 

I      "          1 

0.937                      1/-187 

O.  65~4 

Z>    ... 

•                i    .........  . 

0.999     i    0.6.93     I/.Z6Z     1   0.6,54 

1 

;  I.-2.6.1 

tf.  999 

C-6,74. 

j 

I    ' 

E-Z        ■   o.8U 

■ 

0.06,8     '  (c.934)  1 

1 

/  Zo 

&.(*  70 

\ 

F-4-          \e>.  858 

0.06,8 

0.?i6     [    0.  6,09 

1  fl.  «60 
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ANALYSIS  OF  EXPERIMENTAL  DATA 

Strain  Data  Analysis 

The  moduli  of  elasticity  were  calculated  from  the  strain-pressure 
ratios  using  the  formulas  presented  in  Table  8.  The  theoretical  basis 
for  these  formulas  is  shown  in  Appendix  C.  A  summary  for  the  calcula- 
tions for  E^  and  E^   is  given  in  Tables  9  and  10.  No  values  were  cal- 
culated for  the  B-model  because  premature  failure  of  the  instrumented 
model  in  this  series  prevented  the  taking  of  sufficient  strain  data. 
The  values  of  E^  and  E^  used  for  the  B-raodels  are  averages  of  the  values 
for  the  other  models. 

Since  available  instability  formulas  apply  to  isotropic  materials, 
it  was  necessary  to  develop  a  method  for  converting  the  principal  mod- 
uli ftp  and  ^  into  an  effective  modulus  EQ.  This  was  accomplished  by 
first  deriving  an  expression  for  the  shear  modulus  of  an  orthotropic 
material.  The  derivation  followed  a  method  used  by  Timoshenko^  (pp. 
54-57)  in  deriving  the  expression  for  the  isotropic  case  and  is  presen- 
ted in  Appendix  D.  The  shear  modulus  thus  found  is 

G  =  ExE^ 

(Ex  +  E,)(l  +*0 
Note  that  for  E^  =  E^  this  expression  reduces  to  the  isotropic  case 

0=    E 
2(1  4^) 

In  the  derivation,  Poisson's  ratio  was  assumed  to  be  constant.  As  noted 

later,  this  assumption  does  not  lead  to  an  appreciable  error.  Then,  to 

obtain  an  effective  modulus,  the  orthotropic  expression  for  G  was  equated 

to  that  for  the  isotropic  case  as  follows: 
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Em 
Let  Ex  =  IcEm  ,   then     G  =  - .        Equating  this  to 


G  =    ^       gives 
2(1  +  */) 


Ee  -  2E^ 


1  +  k 
The  calculated  effective  values  for  the  modulus  are  summarized  in  Table  11. 

In  addition,  the  authors  undertook  the  derivation  of  an  expression 
for  the  buckling  pressure  of  an  orthotropic  cylinder.  This  derivation  is 
shown  in  Appendix  E,  and  the  buckling  pressure  is  given  in  the  form  of  a 
determinant.  The  calculation  of  buckling  pressure  using  this  determinant 
is  very  cumbersome .  Since  the  expression  was  derived  shortly  before  the 
completion  of  the  thesis,  no  evaluation  of  the  determinant  was  made  with 
respect  to  the  experimental  results  contained  herein.  It  is  believed, 
however,  that  it  may  have  merit  in  computing  buckling  pressures  for  ortho- 
tropic  cylinders.  Also,  it  might  be  used  to  verify  the  assumption  of  an 
effective  modulus  for  use  with  the  isotropic  equations.  It  is  noted  that 
the  determinant  derived  for  the  orthotropic  case  reduces  to  the  expression 
for  the  infinite  tube,  with  E^  the  only  modulus,  when  L  is  made  to  approach 
infinity.  This  suggests  that,  in  reality,  an  effective  modulus  for  use  in 
the  buckling  equation  should  depend  on  the  length  of  the  cylinder.  It  is 

also  noted  that  the  determinant  reduces  to  that  given  by  Timoshenko^ 

E 
(p.  449)  when  &%   is  made  equal  to  E^  and  G  is  put  equal  to  , 

2(1  +is) 

except  that  the  term  -  — —  is  added  to  the  last  term  in  row  3, 

h    2 

column  C  to  take  care  of  the  end  thrust. 

Details  of  the  assumptions  employed  in  the  calculation  of  an  effec- 
tive modulus  follow. 

(1)  As  previously  stated,  Poieson's  ratio  was  assumed  to  be  constant. 
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On  the  basis  of  data  provided  by  the  David  Taylor  Model  Basin 
and  the  Zenith  Plastics  Company,  this  value  was  assumed  to  be 
0.15.  The  effect  of  the  variation  in  Poisson's  ratio  on  the 
value  of  the  calculated  modulus  is  small.  For  example,  using 
V  =  0.10  instead  of  0.15  in  the  calculation  for  E^  raises  the 
calculated  value  about  3%   for  both  the  A  models  and  the  B  models. 

(2)  The  stress  field  in  the  model  was  assumed  to  be  in  accordance 
with  the  Lame  thick  shell  theory,  and  the  strain-pressure 
curves  were  assumed  to  be  linear  up  to  the  point  of  collapse. 
As  a  consequence  of  this  assumption  the  bending  stress  due  to 
the  simply  supported  tube  of  finite  length  is  neglected.  As 
will  be  shown  later,  the  bending  stress  prior  to  buckling  is 
a  very  small  quantity. 

(3)  The  effect  of  radial  strain  was  neglected  in  calculating  E^z 
and  Ex.  This  effect  is  small  as  is  shown  below: 

The  radial  strain  on  the  outer  surface  is  given  by:  (see  Appendix  c) 


<«*  =  P 


r_J z//K,/z     KxV 

Br,  ^~TT  E*  ) 


Thus, 


ho„  = 


_   _!_  -j,(]<l/*-  +      K 


It  is  seen  that  the  greatest  error  would  occur  for  the  largest 
value  of  «rz »  which  value  would  correspond  to  the  thickest 
cylinder,  i.e.,  the  cylinder  with  the  least  value  of  K-\,     For 
the  E-model,  Kx  =  6.5;  K2  =  Kx  -  1  =  5.5;  E^  =  5.5  x  106; 
Ex  =  U   x  10°;  V-   0.15;   then  assuming  a  nominal  value  for 
Ej.  =  3  x  106,  m^=  1/3  -  (0.15)(2a2£  +  liSJ  S  0.0^6 
Now  the  exact  formula  for  computing  E  using  outside  strain 
gage  data  is-.      _    _Kt_0  +  f)('~Zl') 


'0    ^ttii  -v)+-2/nnXx  +  z/ma. 


It  was  found  that  the  difference  between  the  value  of  Ey 
calculated  by  using  this  equation  with  n_  =  0.04.6  and 
that  calculated  by  neglecting  a^.  entirely  was  less  than  1%. 
It  was  further  noted  that  the  longitudinal  compressive  strain  was 
greater  on  the  outside  than  on  the  inside  of  the  models.  This  difference 
is  attributed  to  the  fact  that  some  bending  of  the  shell  undoubtedly 
occurs  even  for  the  lower  pressures.  That  this  difference  is  not  inherent 
in  the  strain  relations  can  be  shown  by  the  following  development: 
From  Appendix  C, 

6X2  =  f~  -  v(?^-+  -£-)    (radial  stress  =  r  ) 
6X|=  2 Ex  —  ^Vr-  ' J  (radial  stress  =  0  ) 

Since  K2  =  K^-l, 
e**      t^~2£K  E*      ^~EV        Er  2£"x     +     £p 


€xz-*X,  = 


from  which, 


>^>v,  -  mnY.    = 


/- 


'XX 


x,    =  V{£0 


-to 


Since   E^  *   5  x  106    E,  =3x  106     and  J>  =  0.15, 
mX2  -  mX(  S  (0.15) (1/5  -  1/3)  =  -0.02 

This  difference  is  small  in  comparison  with  the  differences  actually 
measured  and  is  of  opposite  sign.  It  is  concluded,  therefore,  that  the 
measured  difference  results  from  bending  of  the  generators.  In  order  to 
reduce  the  errors  arising  from  this  effect,  an  average  of  mXl  and  mXz 
was  used  in  all  calculations  which  called  for  either  value. 
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Table  7  shows  that  even  for  identical  models  there  were  considerable 
differences  in  measured  values  of  strain-pressure  ratios.  These  differences 
may  be  attributed  to  two  causes: 

(1)  The  properties  of  the  individual  models  may  vary  slightly,  even 
though  cut  from  the  same  tube. 

(2)  Even  on  the  same  model,  the  strain-pressure  ratios  will  vary  from 
point  to  point  around  the  circumference,  as  shown  in  reference  (38). 

The  most  accurate  data  for  use  in  computing  moduli  would  be  that  taken 
from  a  number  of  strain  gages  spaced  around  the  circumference  of  each 
model,  using  the  average  values  of  strain-pressure  ratios.  Due  to  limi- 
tations discussed  previously,  this  procedure  could  not  be  followed. 
Instead,  the  averaging  process  was  accomplished  by  using  inside  and  outside 
strain  data,  and  averaging  the  moduli  thus  calculated.  Where  circumferential 
variations  in  strain-pressure  ratios  are  caused  by  the  tendency  to  lobe 
formation,  the  use  of  inside  and  outside  data  from  w back-to-back"  gages 
tends  to  cancel  the  variation.  In  the  case  of  the  A,  D,  and  E  geometries, 
two  instrumented  models  of  each  were  tested  in  order  to  provide  more 
data  for  averaging. 

Since  the  tubes  were  fabricated  using  the  same  techniques  and  material, 
one  would  expect  the  moduli  to  be  about  the  same  for  all  the  models.  Such 
is  the  case,  except  that  the  A  models  were  found  to  have  a  slightly  higher 
modulus. (See  Tables  10  and  11).  As  previously  noted,  the  effective  modulus 
for  the  B  models  was  taken  as  the  average  of  that  of  the  other  models,  due 
to  insufficient  strain  data.  The  small  differences  in  resin  content  for 
the  models  (Table  11)  are  not  sufficient  to  account  for  the  differences 
in  ~oduli. 
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TABLE  8 
SUMMARY  OF  THICK  SHELL  FORMULAS  FOR  CALCULATING 
MODULI  OF  ELASTICITY 


Inside  gages,  dummy  block  not  subject  to  pressure: 

L?  =  z L*  - 


***,+■  Wx,  Z.(myt-hZlYr\pt) 


Outside  gages,  dummy  block  not  subject  to  pressure: 

t-t- ~, ^ / ^        t.x- 


'n*t(i-v)-hv(yyixi+YyiY)  '  2 [ho X2 (■-*)+ ^7*^  +  171,)] 


Outside  gages,  dummy  block  subject  to  hydrostatic  pressure**: 

tf  =   -  tLx—       ;      ; 

^^i4-^»?xz.  Wis.  +  vwti. 


**  Dummy  block  assumed  to  be  of  same  material  as  model;  circumferential 
active  gages  balanced  against  circumferential  dummy  and  longitudinal 
active  gages  balanced  against  longitudinal  dummy. 
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TABLE  8  (Continued) 
SUMMARY  OF  THICK  SHELL  FORMULAS  FOR  CALCULATING 
MODULI  OF  ELASTICITY 


E$  =  Circumferential  modulus  of  elasticity 

£x  =  Longitudinal  modulus  of  elasticity 

R  =  Mean  radius,  f/z(R,  +  Ra) 

-A  =  Wall  thickness,  (R^-R,) 


Ri-Rf- 

Rl-Rf 

m#,    =  Avg.  strain-pressure  ratio,  inside  circumferential  gages, 
dummy  not  under  pressure. 

rr)x,    =  Avg.  strain-pressure  ratio,  inside  longitudinal  gages, 
dummy  not  under  pressure. 

A7<*2  =  Avg.  strain-pressure  ratio,  outside  circumferential  gages, 
dummy  not  under  pressure. 

/n*2  =  Avg.  strain-pressure  ratio,  outside  longitudinal  gages, 
dummy  not  under  pressure. 

fv)\r    =  Avg.  strain-pressure  ratio,  radial,  on  outer  surface. 

m^t   =  Avg.  strain-pressure  ratio,  outside  circumferential  gages, 
dummy  under  hydrostatic  pressure. 

hOxi   =  Avg.  strain-pressure  ratio,  outside  longitudinal  gages, 
dummy  under  hydrostatic  pressure. 

2^   =  Poisson's  ratio. 
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table:    9 


DIMENSIONS       OF     TU3ES      FROM     Wh/ICH       MODELS       Wf/<fc"      TAKEN 


akID      CALCU  LAT  ION      FOR        K, 


TU3E 

A 

6 

c 

D 

£ 

A\S  INSIDE  D/A  (IN.) 

3.1  fS 

3.  US 

3.  1  ?-T 

3   IIS 

3   I9S- 

-A 

AVG.  WALL  THICKNESS  (IN) 

O.  l<o9 

o.Zos- 

0Z4O 

OZ8I 

0-3  09 

D 
AVGr.    MEAN  VIA.  (/N ) 

3.3t>4 

J.40O 

3.  43S 

3.476 

3S04 

-A/ 

/D                          0.0SoZ4 

0.0  6  02.  J 

0.oi>987 

C08O84 

0. 088 18 

.AVS.  INSIDE  RADIUS  (IN.) 

| 

I.S97S 

I.S97S 

l.S?7S 

/.  Sf  7S 

/.S97S 

1 

AVG.  OUTSIDE  RAD/US  (in.) 

1.  766-S" 

/■  BCZS 

IS37S 

1.876* 

I.706S- 

7.  1Z0S 

3.Z4  7  0 

3-37i>4 

3.  SZ88 

3.0,347 

*: 

2-  SSI  o 

2   SSZO 

z .sszo 

z-sszo 

2SSZO 

Rl  -  *,' 

O.S(,8S~ 

O  6.970 

0.8Z.4  4 

0-  97  G  8 
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Stress  Analysis 

In  order  to  properly  evaluate  experimental  results,  it  was  necessary 
to  verify  the  mode  of  failure  for  each  model.  Although  all  models  were 
designed  to  fail  by  instability,  the  predicted  instability  failure  pressure 
for  the  D  and  E  models  was  fairly  close  to  the  corresponding  compressive 
yield  failure  pressure  predicted  using  the  thick-shell  theory  and  a  nom- 
inal value  for  the  yield  strength  of  the  material.  Also,  the  thick-shell 
formulas  used  to  calculate  the  stress  do  not  take  into  account  the  bend- 
ing effects.  Because  of  these  considerations,  the  following  steps  were 
taken  in  order  to  make  an  evaluation  for  the  mode  of  failure. 

(1)  In  order  to  evaluate  the  effects  of  bending,  a  formula  was  de- 
rived to  calculate  the  longitudinal  bending  stress  for  a  simply 
supported  cylinder.  The  method  of  derivation  was  suggested  by 
Mr.  John  Pulos  of  DTMB  and  is  presented  in  Appendix  F.  The  ex- 
pression for  the  bending  stress  at  midbay  is: 

Substituting  values  for  the  E  model,  the  theoretical  longitud- 
inal bending  stress  at  midbay  was  found  to  be  only  75  psi. 

(2)  The  Lame  (thick-shell)  stress  corresponding  to  the  failure  pres- 
sure was  calculated  for  each  of  the  models  and  compared  both 

with  a  nominal  compressive  strength *of  70,000  psi  as  obtained 

33 

from  DTMB  Report  LU-3  ,  and  with  the  lowest  calculated  value 

obtained  from  a  yield  model  (63,233  psi).  For  the  E  model,  the 
highest  calculated  stress  was  59,285  psi  (see  Table  6).  This 
value  is  10,715  psi  below  the  value  for  nominal  strength  but 
only  3,948  psi  below  that  obtained  from  the  yield  model.  It 

*  Yield  strength  in  circumferential  direction,  o^ 
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should  be  noted,  however,  that  application  of  the  Hencky-von 
Mises  failure  criterion  would  cause  the  differences  to  be  greater. 
(3)  Finally,  the  physical  appearance  of  the  collapsed  models  was 
compared  to  that  of  the  thinnest  model,  for  which  a  pure  in- 
stability failure  was  certain,  and  with  that  of  a  short  model 
which  definitely  failed  by  yielding  of  the  material.  Figures  11 
and  12  show  collapsed  C,  D,  and  E  models  together  with  a  short 
model  which  fa: led  by  yielding.  It  is  noted  that  the  longer 
models  pictured  failed  along  a  longitudinal  line  parallel  to 
the  model  axis  by  buckling  inward,  whereas  the  short  model  tended 
toward  failure  throughout  most  of  its  wall  area. 
On  the  basis  of  these  considerations  it  was  determined  that  the  A, 
B,  and  C  models  failed  purely  in  the  instability  mode.  Since  the  com- 
pressive stress  for  the  E  models  was  within  about  4000  psi  of  the  value 
obtained  from  a  yield  model,  there  still  may  be  reason  to  question  whether 
or  not  some  yielding  took  place.  Also,  there  was  scatter  in  the  points 
for  the  D  and  E  models  (See  Fig.  13)  which  may  be  a  result  of  local  yielding 
caused  by  stress  concentrations  due  to  the  heterogeneous  nature  of  the 
material. 
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PRESENTATION  AND  nTSm^Tpg  0F  BMTtTTg 

Result*  of  the  exponents  are  shown  graphically  in  Flgure  „  The 
critical  pressures  deterged  f™  tests  were  converts  tc  the  „on-dlBeneional 
'o™  P/Se  .  These  value8  „,  plotted  agalngt  ^  thlckness_diameter 

ratio.  Calculations  for  p/Be  are  shown  in  Table  12. 

Also  plotted  on  Figure  13  are  curves  representing  various  theoretical 
buckling  formulas,  using  an  average  value  of  L/D  for  the  models  tested. 
(Average  L/D  =  2.53). 

in  spite  of  the  anisotropic  nature  of  the  n,odel  material  and  the 
inaccuracies  inherent  in  the  computation  for  an  effective  modulus,  very 
good  agreement  with  theory  was  obtained  e^for  the  thicker  models.  The 

that  the  thin  shell  buckling  theory  is  not  adequate  for  ratios  of  h/D  i„ 
this  range.  As  previously  noted,  however,  so„e  yielding  My  have  occurred 
in  these  models. 

It  appears  that  the  von  Mi...  equation  is  the  best  one  to  use  in  the 
-nge  of  h/D>  0.05  since,  in  general,  it  yields  conservative  values  for 
the  collapse  pressure.  *sig„  curves  for  the  buckling  of  si.ply  supported 
cylinders  were  therefore  drawn  up  for  the  von  Mises  equation.  (See  Pig.  U) 
Data  for  these  curve,  were  calculated  using  the  IBM  7090  counter  at  DTMB. 
These  data  are  reproduced  in  Appendix  0. 

Since  curve,  were  not  available  for  deterffli„i„g  the  number  of  de- 
ferential lobes  into  which  a  relatively  thick  cylinder  would  fail,  curves 
were  drawn.   (See  Fig.  15).  These  _„  „„  ^^  ^  ^  ^ 

Appendix  G  by  putting  the  buckling  equation  in  the  for* 

P/K  =  Fjfh/D)3  +  F2(h/D) 
for  each  n  and  L/D  and  then  solving  for  the  desired  intersections.  These 
curves  are  for  z;  =  0.I5. 
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For  values  of  n  with  V  =  0.3»  in  the  range  h/D  £  0.02,  aee  reference 

(33). 

As  far  as  is  known  by  the  authors,  these  curves  for  the  von  Mises 
equation  and  for  the  number  of  lobes  are  the  only  ones  available  in  the 
range  of  h/D  greater  than  0.05. 
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CONCLUSIOHS 

Based  on  the  results  of  these  experiments  the  following  conclusions 
are  drawn: 

(1)  The  von  Mises  equation,  equation  (6),  Ref.  22  is  the  best  in- 
stability equation  for  cylinders  with  ratios  of  h/D  greater  than 
0.05.  It  gives  conservative  values  for  the  buckling  pressure  up 
to  an  h/D  of  about  0.07.  Use  of  presently  available  simplified 
forms  of  buckling  equations  cannot  be  expected  to  give  conserva- 
tive values  for  the  buckling  pressure  in  this  range  of  h/D.   It 
is  noted  that  Windenburg  and  Trilling^  had  previously  reached 
the  conclusion  that  the  von  Mises  formula  was  probably  the  best 
one;  however,  their  conclusion  was  based  on  an  analysis  of  exper- 
imental and  theoretical  results  for  cylinders  with  ratios  of  h/D 
less  than  0.007  in  contrast  with  the  much  thicker  range  of  0.05 
to  0.088  investigated  in  this  thesis. 

(2)  For  practical  cases,  the  use  of  stiffening  rings,  heavy  webs,  and 
bulkheads  would  result  in  smaller  ratios  of  h/D  and  L/D  for  a 

given  weight  of  structural  material.  It  is  concluded,  then,  that  for 
most  practical  cases,  the  thin  shell  theory  is  sufficiently 
accurate  for  predicting  instability  failures. 

(3)  For  cylinders  built  up  of  fiber  glass  reinforced  plastic  layers, 
the  methods  employed  in  this  thesis  are  sufficiently  accurate 

to  predict  instability  failure  pressures. 
(4.)  For  an  orthotropic  cylinder,  the  effective  modulus,  as  used  in 
the  buckling  equation,  is  probably  a  function  of  the  length 
between  supports. 
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RECOMMENDATIONS 

Further  Investigations, 

The  experimental  apparatus  designed  and  assembled  for  this  thesis 
provides  a  means  for  the  undertaking  of  a  variety  of  experimental  in- 
vestigations in  the  field  of  hydro statically  loaded  shells.  Some  investi- 
gations which  might  be  undertaken  are: 

(1)  Investigate  the  stresses  in  short,  relatively  thick,  fiber- 
glass reinforced  cylindrical  shells.  Measure  the  strain 
distribution  along  a  generator  and  determine  the  effects  of 
bending.  Verify  the  formula  for  bending  stress  which  is 
derived  in  Appendix  F. 

(2)  Using  fiber-glass  reinforced  plastic  models,  verify  the  ortho- 
tropic  buckling  determinant  as  derived  in  Appendix  E.  In 
conjunction  with  this,  evaluate  the  determinant  to  arrive  at 
a  general  expression  for  buckling  pressure  similar  to  the 
von  Mises  equation  for  the  isotropic  case.  Using  the  result, 
develop  a  method  of  computing  an  "effective  modulus"  for  use  in 
isotropic  buckling  formulas. 

(3)  Perform  experiments  on  models  of  unusual  form,  such  as: 

a.  Stiffened  spheres 

b.  Ellipsoids 

c.  Sphere-cylinder  intersections 
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RECOMMENDATIONS  (Continued) 

improvement  of  Apparatus, 

Replace  valve  (G) ,  (Fig.  U)   with  one  rated  at  10,000  psi  or  higher, 
such  as  American  Instrument  Company  valve,  Cat,  No.  44.-1505. 

Install  a  fitting  in  the  line  to  each  gage  to  damp  out  the  sudden 
drop  in  pressure  when  a  model  breaks,  thus  protecting  the  gage  pointer. 
There  are  commercial  fittings  available  for  this  purpose. 

If  a  greater  number  of  strain  gage  leads  are  required  to  be  run  for 
the  gages  mounted  on  the  outside  of  the  model,  it  is  recommended  that  a 
CONAX  Corporation  thermocouple  fitting,  Cat.  No.  TG-20-B8,  be  installed 
in  place  of  the  fitting  "A"  nov  installed.  The  TG-20-B8  will  take  eight 
#20  wires. 
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APPENDIX  A 
RECORD  OF  STRAIN  GAGE  DATA 

The  following  tables  consist  of  the  recorded  strains  corresponding 
to  the  gage  pressures  and  corrected  pressures.  The  model  number,  gage  type, 
gage  factor,  date  of  test,  gage  location  ,  cycle  number,  run  number,  and 
pertinent  remarks  are  also  indicated. 

Three  types  of  Baldwin-Lima-Hamilton  gages  were  used  for  the  strain 
measurements: 

Gage  Type  Nominal  Gage  Length 

A-5-1  1/2" 

FAP-50-12  1/2" 

FA-100-12  1" 

The  type  A-5-1  gage  was  used  on  the  inside  of  the  model  only,  since 
the  wire  gage  is  particularly  sensitive  to  hydrostatic  pressure.  The 
FAP-50-12,  a  foil  gage,  was  used  both  inside  and  outside  the  models.  The 
FA-100-12  was  used  on  the  outside  of  model  k-U   only.  All  three  types 
appeared  to  give  reliable  data.  No  particular  difficulty  was  experienced 
in  handling  and  applying  the  gages. 

In  the  row  entitled  Gage  Location.  "L"  stands  for  a  longitudinally 
mounted  gage,  nC"  for  a  circumferentially  mounted  gage,  "1"  for  an  inside 
gage,  and  "2"  for  an  outside  gage. 

In  some  cases  the  gages  were  read  one  at  a  time  requiring  a  differ- 
ent pressure  cycle  for  each  gage.  A  record  was  kept  of  the  total  number 
of  cycles  for  each  model  and  the  number  of  runs  for  a  particular  gage. 
These  data  were  taken  in  order  to  ascertain  the  effects  of  creep,  if  any. 
It  appeared  from  the  results,  however,  that  creep  was  not  an  important 
factor  for  these  short  terra  tests. 
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APPENDIX  B 
PLOTS  OF  PRESSURE  VERSUS  STRAIN 

This  section  contains  plots  of  the  measured  strains.  The  resulting 
slopes,  designated  by  the  letter  n  with  appropriate  subscripts  are  indi- 
cated on  the  plots. 

It  is  noted  that  all  plots  are  very  nearly  linear  until  the  col- 
lapse pressure  is  approached.  The  reversed  slope  for  gage  L-2,  model 
A-4  indicates  the  development  of  a  tensile  stress  as  a  result  of  bend- 
ing near  the  collapse  pressure. 
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APPENDIX  C 

THEORY  FOR  CALCULATION  OF  E^  AND  Ex  FROM  STRAIN  DATA 

The  circumferential  stress  in  a  thick-walled  cylinder  under  external 
hydrostatic  pressure,  p,  is  given  by: 

to  rue  point 

where  r  is  the  distance  from  the  centerAat  which  the  stress  is  measured. 
If  r  is  given  values  of  r  =  R,  and  r  =  R2,  the  following  expressions 
evolve  for  the  circumferential  stress  on  the  inner  and  outer  walls  of 
the  cylinder: 

(<j^  and  O^are  positive  if  compressive.) 
The  longitudinal  stress  is  constant  over  the  cross  section,  and 
is  given  by: 

<rx-  rj^r  (tot.  42) 

Letting  K   -  1      we  have 

Rt-R,1' 
<r*  =   rK'M  ) 

If  we  assume  a  material  which  has  three  principal  values  of  Young's 
Modulus  and  a  single  value  of  z^,  then  from  Hooke's  Law: 
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_  G± z^OV  _  ±>  CTX 

^0      u      C-x 

£*      Cf      £,!5 

Combining  Hooke's  Law  and  thick  shell  theory,  ve  have  strains  in  the 
inside  and  outside  surfaces  of  a  thick-walled  cylinder  subject  to  external 
hydrostatic  pressure: 

0i   *VE*      E,      ZE* 

'K,     7J  vK-i. 


t 


X2. 


(Note  that  0;=  P  on  the  outside  surface) 

**•  =r(it;~itj  (we)      I     (4.) 

(Mote  that  CV=  0  on  the  inside  surface) 

The  stresses  in  the  surfaces  of  a  block  of  material  subject  to 
external  pressure,  p,  are  given  by: 

°V  —  <Tx  =  0\-  -  ~p 

From  Hooke's  Law,  the  strains  on  the  outside  surface  of  the  block 

are  given  by: 

/  /   _  _j-j_ & 

\t^\rD  t,XJ>  Efij> 

The  d  subscript  on  a  term  indicates  that  it  refers  to  the  block  rather 
than  to  the  cylinder. 
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Consider  the  case  of  strain  measurements  made  with  strain  gages 
on  the  inside  surface  of  the  cylinder,  using  a  dummy  gage  not  subject  to 
pressure.  The  strains  as  read  are  given  by: 

Solving  the  first  two  of  these  equations: 
K,(i  -vx) 


T)£  -f-Z>  tY)X, 

K,  (i  -v2) 


to 


z{mnx^vrY\^ 


Where      ^   -  6*./y? 


Therefore  a  model  instrumented  thus  can  be  subjected  to  hydrostatic 
pressure  and  the  measured  strain-pressure  ratios,  m  and  mpf  ,  with 
an  assumed  or  known  value  of  is   can  be  substituted  to  solve  for  E^.  and  E^, 
The  radial  strain  does  not  enter  the  solution,  and  the  radial  stress  is 
zero. 

Consider  the  case  of  strain  measurements  made  with  strain  gages 

on  the  outside  surface  of  the  cylinder,  using  a  dummy  gage  not  subject 

r 

to  pressure.  The  strains  as  read  are  given  by: 

/Kz  _  v      _     isK, 

Solving  this  system  of  equations  simultaneously  gives: 
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E*   = 


Ex  - 


(9) 


^[^xz(i-^+-  z^m^-^rj] 


Consider  the  case  of  strain  measurements  made  with  strain  gages 
on  the  outside  surface  of  the  cylinder,  using  dummy  gages  mounted  on  a 
block  subject  to  the  hydrostatic  pressure,  the  longitudinal  active  gage 
being  balanced  against  a  longitudinally  oriented  dummy,  and  the  circum- 
ferential active  gage  being  balanced  against  a  circumferentially  oriented 
dummy.  Assume  the  dummy  block  to  be  of  the  same  material  as  the  cylinder. 

The  indicated  strain  readings  are: 

where  the  primes  indicate  that  the  strains  are  the  indicated  values  rather 
than  the  actual  strains.  From  (3)  and  (5): 
'  K  2.    tj         *sK, 
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Eis 


(fO) 


6^      =      0 


C* 


The  terms  involving  radial  strain  drop  out.  A  comparison  of  (10) 
with  (3)  shows  that  they  are  identical,  except  that  in  (10)  the  radial 
termd  are  missing,  and  K2  -  1  and  K]/2  -  1  appear  in  place  of  K2  and 
Kn/2,  respectively.  Thus  by  having  the  compensating  gages  arranged  as 
in  this  case,  the  values  of  E^  and  E  may  be  computed  from  the  indicated 
strains  6^  and  <=Xl,  and  the  radial  strains  do  not  enter  the  solution. 

Solving  (10): 


£*  = 


^a.  +  ^  m'X3L 


00 
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APPENDIX  D 
DERIVATION  OF  SHEAR  MODULUS  FOR  AN  ORTHOTROPIC  MATERIAL 

Following  a  method  similar  to  that  used  by  Timoshenko  (Ref.  UO, 
pp.  54-57)  for  developing  the  shear  modulus  for  the  isotripic  case,  a 
formula  is  derived  for  the  shear  modulus  of  an  orthotropic  material 
by  assuming  two  principal  moduli  of  elasticity  and  a  constant  value 
for  Poisson's  ratio. 

Consider  an  element  in  pure  shear.  For  this  condition  07  =  -  *&  =  ?Ze 
as  is  readily  apparent  from  Mohr's  circle.  The  element  can  be  represented 
as  shown  below: 


<re 


<rK 


1  i  1  1  I  M- 


61 


r  i  n  t  n  t 


<ra 


(The  deformed  element  is  represented  by  dotted  lin^s.)  The  deformed 
angle  at  A  is^r-V;  that  at  B  is  -f  +  y  ,  where  Y  is  the  shearing  strain. 

A/our       ca'  —  OA  ( 1+  e„)        *»</        &s'  =   o~3  it*  €*) 


£«   - 


^C7S 


6  --S 


(Tk 


-Dl  - 


hut         <TK   ~     TC&  d"^  <7©    =   _    ? 


<& 


7~h  <=  >? 


0/1 


—     "  i  'r"&         f  r*»~\  its'        ttx\  i       2**  t.  ^l 


S/nce.      OA    =  03 


_[/- 


^6  __    ^^, 


3 


M/v  %.  -  ta*  7z  __    /  -  Vz 


SoJv/nj      for       Y~    7 


/  VCS 


Y  = 


AS£±lI±KlU!± 


2E«E*    __  (EK-Ee%l-^) 


A/e<j /ecT/ h<j      7U<s     Second     fe.r-*%»     rn     the    cfcnom/  axp  To  r-      Ahd     So/r/r*<j      for     7a^4 
Ek  C  g ^       |^ 

Then        the.     o*-thotrop/<-     c  x.  press /o*>      far-      G-     become  s 


G-  = 


Al 


E«  Eg, 


(E<  +  E0)(/^-^J 


Conway   arrived  at  a  similar  expression  using  an  energy  analysis. 
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APPENDIX  E 

DERIVATION  OF  EQUATION  FOR  THE  CRITICAL  PRESSURE 
OF  A  SIMPLY  SUPPORTED,  ORTHOTROPIC,  CIRCULAR 
CYLINDER  SUBJECTED  TO  HYDROSTATIC  PRESSURE 

Stresa-strain  relations. 

For  isotropic  shells,  the  relation  between  the  stress  resultants 
and  couples  and  the  corresponding  strains  is: 

MK  =  TT^T  (€*o   ^^eo)  MK   =  ~D  (Xx   *  ^Xe) 

N^^T^-  (6eo  +  *€*-)  M6  =  -D(X0  +  *Xk) 


where  the  subscript  "0" refers  to  strains  in  the  middle  surface. 

If  is  is  assumed  that  the  cylinder  is  orthotropic  with  principal 
moduli  Ex  and  Ee  and  with  a  constant  value  of  Poisson's  ratio,  these 
equations  may  be  written: 

M*  =  JT^I  (*"•€#•* *£e,0J  Mi  =  -  jtcT^)  (£*  X*  t  *E*  z«) 


-El- 


The  strains  and  changes  of  curvature  in  terms  of  the  displacements 


are: 


Ju 

Jx 

CL       d© 

CL 

too   = 


Vo   —  \  v  "*"  cl      ^a 


^>ce  ~"  a.  \   J  x    e)«  J©/ 


The  stress  resultants  and  couples  can  then  be  written  in  terms  of 
the  displacements: 

h 


N*  = 


/-y* 


{a-    6Q  ~cZ) 


/Ve  = 


e         /  -  ?s 


a.    e^©         a, 


-Or   *-  ^^x 


^/~     ^  ^ 


A,' 


M9  = 


~/z(/-ry[_   a-*    U©         <^©*/    *  ^^x    <>** 


Mr.   =    £>> 


xe 


«^X  ^       c)©J 


/V\„«  =  <^ 


'X© 


£*- 


<9x  <^x:  £e 
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The  Equilibrium  Equations . 

The  equilibrium  equations  for  a  circular  cylinder  under  radial 
load  only  are: 


a-    -j-r  +  &   tt  =  0 


d  x 


c)  0 


d  X 


e)Al 


»© 


ax 


=  o 


=    O 


See  Timoshenko,  Ref.  30,  p. 447;  note  that  the  term  yo_  \~J7Je        d  x.  J 
has  been  neglected,  as  was  done  by  von  Mises  .  J 

These  equations  may  be  modified  to  include  the  effect  of  axial 
load  as  follows.  This  procedure  was  used  by  von  Mises. 

The  total  axial  force  on  the  ends  of  the  cylinder  is  equal  to  77" CL  p 
This  force  is  transmitted  to  the  ends  of  the  cylinder  of  area  2  rra.fi  . 
Hence  the  pressure  in  the  cylinder  walls  is 

P    —  Zrrc^h  Zh 

How  consider  an  element  of  the  shell  in  the  longitudinal  plane. 


-£3- 


After  deformation,  the  curvature  of  the  generatrix  (neglecting 
higher  order  terms)  is 


The  force  per  unit  length,  p'h,  thus  has  a  component  in  the  a-direction 


of 


2. 


2-       J  x* 


Since  the  radius,  a,  had  been  divided  out  of  the  original  equilibrium 
equation,  this  term  is  written  as 


«-  ^z 


p  (  or  y±r) 


In  order  to  include  the  effect  of  the  axial  load,  this  term  is  added 
to  the  term  in  parenthesis  in  equation  (3)  to  gett 


Boundary  conditions. 

The  transverse  and  longitudinal  cross  sections  of  the  cylinder  assume 
the  following  shapes  after  buckling: 


—  X 


h~  ^ 


*»/ 
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At  x  =  0,  the  displacements  may  be  written: 
U   =  A  s/n  n  e     \ 

/V  —  B  cos  n  e     )      a/"  /s     90°  out    of   phase.    <*<th    u    at  J  -W. 
<cw  -  C  sin   ai  ©  / 

At  x  =  0,  u=0  and  v  and  w  are  maximum;  hence,  the  displacements  for 
any  value  of  x  become: 

u   =  A  s*t*  n  ©  s/"  ~J~ 

/\Z"  —   O   cos  r>Q   cos    ~jT 
ajT  -   C  s/n  hB   cos  -jf 

Let         \  =    "jf"    ,     then 

U  =  A  s/n  n  o  s/a/  /a.* 
/vr  =  B  cos  »e  cos   A.  K 

AaT  -  C  J//V  A»0  cos    /a.  * 


Determinant  for  the  Critical  Pressure, 

By  substituting  the  boundary  conditions  into  the  relations  for  the 
stress  resultants  and  couples,  and  substituting  the  resulting  expressions 
into  the  equilibrium  equations,  a  system  of  three  equations  in  A,  B,  and 
C  is  determined.  In  order  that  an  instability  condition  exist,  the  constants 
A,  B,  and  C,  which  determine  the  magnitudes  of  the  displacements,  must 
have  values  other  than  zero.  This  condition  requires  that  the  determinant 
of  the  coefficients  of  A,  B,  and  C  be  zero.  The  determinant  for  this  case 
is  given  on  the  next  page.  Evaluation  of  the  determinant  for  a  given 
geometry  and  material  yields  the  critical  buckling  pressure. 
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APPENDIX  F 

DERIVATION  OF  AN  EXPRESSION  FOR  THE  LONGITUDINAL  BENDING  STRESS 

Hi  L 
SIMPLY  SUPPORTED  CIRCULAR  CYLINDER  SUBJECTED  TO  EXTERNAL  HYDROSTATIC 

PRESSURE 

Salerno  and  Pulos43  have  derived  expressions  for  the  stresses  in 

a  circular  cylindrical  shell  supported  by  uniformly  spaced  circular 

rings  of  constant  cross  section.  By  using  the  displacement  function 

obtained  in  their  paper  with  the  appropriate  boundary  conditions  for 

simple  support,  an  expression  for  the  longitudinal  bending  stress  at 

midbay  is  derived. 

The  longitudinal  bending  stress  is  given  by: 

Eh  d^±r 

The  displacement  function, ><v-,  as  derived  by  Salerno  and  Pulos  for 


deflection  in  the  radial  direction,  is 
sur  =  B  cosh  (A,  x)  +  F  cosh  (As  x)  — 


Eh 


where  B  and  F  are  constants  depending  on  the  boundary  conditions,  X,  ¥*  A3)  and 


**  =      [_"  42    *  1K    +o)    ~~\DR*-)  _ 


A,=      - 


PR     At   pry    (Ah  ) 

4D  '   \\+Dj    ~\DRV 


-|& 


Eh 


D    ~  /£  ( I  -  is*)  (  Fl  CKural     rif'J/Ty     of     she//) 


R    =     Me sh    radivs     of    she//. 


-  Fl  - 


Taking  the  x-axis  along  the  axis  of  the  cylinder,  with  the  origin  at 
the  mid-length,  the  boundary  conditions  for  simply-supported  ends  require 
that 


=  Vz 


1^0      i+       *  =  L/z 


^ur  -    O         a.f 
A  X 


from  which, 


PR 


O  =  B  cos*  A,l/z  +  F  cosh  A3  Vz 

o  -  B A,  cosh  A,  /z  +■  F A3  cos*  A3  /z. 


C'-i) 


Solving  simultaneously, 


»■fcFj')fe)[«""-*,] 


Krenzke  and  Short  **  have  presented  a  graphical  method  for  calculating 
the  stresses  obtained  by  the  analysis  of  Salerno  and  Pulos.  In  this 
reference,  different  parameters  were  employed  as  a  simplification. 
These  parameters,  defined  below,  are  used  for  the  remainder  of  this 
derivation. 


Z6    ,  v 

^  1  =  L     (?>+*  fa) 


A3—      L      (ft  ~  -*  ?*-' 


/  *  r~ 
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Using  the  Krenzke  and  Short  parameters, 


A, 


\ * 


A,*-  A 


(  7> +  +  ?*.£ 

4+  f,  J*. 


A»        _    (7,  z±£*2 
A,z-  A  J        "   4*7, 


7' 7' 


a=- 


iikz^j±L 


7'7i 


PR 


'{'-*/*.) 


cosh  (*7,  e-h^^e) 


F  = 


4-^i  y,  y*_ 


i     PR' 


Eh 


('-%) 


Cos#(y,9-A.pJto) 


Th 


e/i, 


>ur  — 


14*7.7*.     J 


n   PR    /      ^/)  ZOx 


(7>-<7<) 


cosh  (^,6-Xpe) 


4+7'?*. 


,1  PR 


£%(>-%)  cwWfi+ifr) 


COSH  (ft,  3  +*.*7z&) 


Eh 


Taking  the  second  derivative  of  this  ejcpression,  simplifying,  and  sub- 
stituting into  the  expression  for  the  longitudinal  bending  stress  gives: 


<5 


/>**(/-%) 


L  (/-*-*; 


'jLT\It£±jflt 

l  /  I    7>  7*- 


X 


J/wfly/j/<ve^  cosh  ^-£^?,  cojH^—fji 


COSH* Of,  cos^o/pi  +  S'HH^Gf,  S/#*6pi 


cosh  erf,  cos  &f*  s/mh  —£-*?,  s/h  -^~"f* 

COSH     9»,  COS*GI7t   +■  S/HHZe»,  S/HZS/?£ 


At      hniJba<i       *    ~°   ,      J"cJ 


1 


X.-O      L 


rn'O-'AJ 


(/-*■*) 


L  )  Iff. 


J/H/i&V,  S/A  6  ft 


cosH^ep,  cos^eoji  t  j~'mh*&>7,  *'**&'?£ 
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APPENDIX  G 
SOLUTION  OF  VON  MISES'  EQUATION  (6) 

This  appendix  contains  solutions  of  the  von  Mises  equation  (6)^2 
as  computed  by  the  IBM  7090  computer  at  DTMB.  Values  of  p/E  were  com- 
puted for  a  range  of  L/D  from  0.50  to  3.00  in  increments  of  0.25,  and 
for  a  range  of  h/D  from  0.01  to  0.10  in  increments  of  0.01.  The  solution 
for  a  particular  h/D  was  started  with  a  value  of  n  =  2.  The  value  of  n 
was  then  incremented  by  one  until  a  minimum  value  for  p/E  was  obtained. 
Columns  3,  4,  and  5  give  the  auxiliary  quantities  PjMt)   and  zyt  for 
use  in  the  von  Mises  equation. 
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